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Fuzzy Bags and Wilson Lines 
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I start with an elementary observation about the pressure in the deconfined phase of a 
SU(3) gauge theory without quarks. This suggests a "fuzzy" bag model for the analogous 
pressure in QCD, with dynamical quarks. I then sketch how the deconfined phase might be 
described using an effective theory of Wilson lines. To leading order in weak coupling, the 
effective electric field appears in a form familiar from the lattice theory of Banks and Ukawa. 



§1. Fuzzy Bags 
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The spectacular success of the heavy ion programs at the SPS and RHIC jus- 
tifies a careful analysis of the phase transition(s) of QCD at nonzero temperature. 
In this Proceeding I summarize some recent work of mine, 1 ) hopefully in a more 
comprehensible fashion. 

Any fundamental understanding of these phase transitions rests upon the bedrock 
provided by numerical simulations of lattice QCD. In this section I begin by looking 
at old data 2 ) in a new way. This was mentioned in a footnote, Ref. 37 of Ref. 
1. While the data is for a pure gauge theory, it shows that except very near T c , 
the critical temperature, there is a exceedingly simple form for the pressure. This 
immediately suggests an ansatz applicable to QCD, and which might be of use for 
phenomenology. 
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Fig. 1. Lattice results ' for the pure SU(3) gauge theory: to the left, (e — 3p)/T ; to the right, the 
same quantity times T 2 . 

Ten years ago, a group at Bielefeld 2 ) computed the thermodynamic properties 
of a SU(3) gauge theory, close to the continuum limit. In equilibrium, while all 
thermodynamic quantities follow from the pressure, it is convenient to plot what is 
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usually called the "interaction measure" , (e — 3p) /T 4 , where e is the energy density, 
p is the pressure, and T the temperature. This is plotted in the left panel in Fig. 1. 
The interaction measure is the trace of the energy momentum tensor, divided by T 4 , 
and so vanishes if the theory is conformally symmetric. A pedestrian way of seeing 
this is to note that e — 3p = T(dp/dT) — 4p; if the theory is conformally invariant, 
the pressure is just a pure number times T 4 , for which e — 3p vanishes. 

Thus the interaction measure is a dimensionless number which quantifies the 
deviation from conformal symmetry. Because the conformal anomaly is proportional 
to the /3-function, this also measures the deviation from ideality. As can be seen from 
the left panel, the interaction measure is very small below T c , rises steeply around 
T c , with a sharp peak at T max ~ 1.1 T c . Above T max it trails off relatively slowly. 
By Tp ert ~ 4.0 T c , its value is, within a factor of two, equal to that expected from 
perturbation theory, 1 ) where the interaction measure begins as ~ a 2 . 

My concern is with the fall off between T max and T port . In the right panel of Fig. 
1, I take (e — 3p)/T 4 , and simply multiply times T 2 . As can be seen, (e — 3p)/T 2 
is essentially constant. For the temperatures shown, then, the pressure is a sum of 
just two pieces: an ideal term, ~ T 4 , and a new, non-ideal term, ~ T 2 . 

If one ignores the overall normalization, the relative normalization of the two 
terms can be computed without further ado. In an asymptotically free theory, at high 
temperature the pressure approaches that of an ideal gas, so at any temperature it is 
natural to compare the pressure to that of the appropriate ideal gas. In the pure glue 
theory, T c ~ 270 MeV is much smaller than the lightest glueball mass, ~ 1.5 GeV, 
and so, relative to the ideal gas, the confined pressure is very small. 2 )' 3 )' 4 ) Thus the 
pressure (nearly) vanishes at T c . This gives 

P P urcgiuc(r)«/ pcrt (T 4 -r 2 r 2 ) , i.it c <t<4.ot c , (i-i) 

where / per t is a constant. 

What of other numbers of colors, N c l The deconfining phase transition is of 
second order 3 ) for N c = 2, and of first order 2 )' 4 ) when N c > 3. As iV c increases, the 
transition becomes more strongly first order, 4 ) with a latent heat ~ N 2 . For three 
colors, 2 ) the transition is weakly first order. 

The change in the order of the transition with N c affects the interaction measure, 
but not dramatically so. 4 ) Since the pressure is continuous at T c , and as the confined 
phase has negligible pressure, then whatever the order, the pressure is almost zero 
at T c , p(T c ) ~ 0. In contrast, the energy is sensitive to the order: while e(T~) ~ 0, 
just above the transition the energy vanishes for a second order transition, and is 
nonzero for a first order transition. 

For a first order transition, then, at T c the interaction measure ~ e(Tj l ")/T 4 . 
For the strongly first order transition with four colors, 4 ) it is not surprising to find 
that the maximum in the interaction measure is at T c = T max , and that it falls off 
after that. 

For a second order transition, the energy is continuous at T c , so e(T c + ) ~ 0, and 
the interaction measure is nearly zero at T c . For the deconfining transition with 
two colors, 3 ) the interaction measure increases from near zero at T c , has a sharp 
maximum at T max ~ 1.15 T c , and then falls off after that. Thus for two colors, 3 ) the 
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peak in the interaction measure can be viewed as a remnant of that for an infinite 
number of colors. 4 ) Since the transition is nearly second order for three colors, 2 ) 
e(T+) is small, and the interaction measure looks like that of two colors; the peak 
in the interaction measure moves closer to T c , to T max ~ 1.1 T r . 

The formula in (1-1) only applies above the maximum in the interaction measure, 
so we should only compare different N c at T > T max . For N c = 4, this is for all 
T > T max = T c . This appears to be supported by lattice simulations. 4 ) To be fair, 
the data for N c = 8 does not, but perhaps lattice discretization errors are larger 
there. 4 ) For two colors, one compares for T > T max = 1.15 T c : the data of Ref. 3 
appears to violate (1-1) by a large amount, ~ 50%. However, it is not clear that 
these simulations are close to the continuum limit, and so new simulations would be 
most welcome. 

Before describing how one might extend (IT) to other temperatures, to T < T max 
and T > T per t, I skip ahead directly to the case of dynamical quarks and three 
colors. For three flavors or less, 5 ) while the "transition" often becomes a crossover, 
an approximate "T c " can still be defined. Whatever the order, though, the interaction 
measure behaves similarly. Relative to an ideal gas of quarks and gluons, there is a 
small but significant pressure below T c , and a maximum in the interaction measure 
at a temperature T max , which is above T c . The surprise about (IT) is that the 
leading correction to the ideal gas term is ~ T 2 , and not ~ T 3 . I speculate that this 
is generic: that for T max < T < T pert , the pressure is a series in powers of 1/T 2 times 
the ideal T term. I call this a "fuzzy" bag model for the pressure: 

PQCD(T) ~ /port T 4 - i?f uzzy T 2 — Byiir + ■ ■ ■ , T max <T< Tp ert (T2) 

The upper bound, T pcr t, denotes the point at which perturbative contributions to 
the interaction measure are of the same order as that found from the lattice; it is a 
few times T c , something like ~ 4T C . 

In (IT), /p er t is dimensionless, Bun is the usual MIT bag constant, 6 ) with 
dimensions of (mass) 4 , and B{ uzzy is a fuzzy bag constant, with dimensions of (mass) 2 . 
For the pure glue theory, from (IT) B {uzzy = / pe rtT 2 , and B M rr < ^fuzzy- 
Recent lattice simulations 5 ) appear to support (T2). Note that with a fuzzy bag 
constant, the sign of the MIT bag constant is not guaranteed, but these simulations 
find that £>mit is positive, 5 ) as in the original MIT bag model. 6 ) 

With a little work, it should be possible to generalize (1-2) to the entire decon- 
fined phase. Above T per t, we can ignore the non-ideal terms, and include only f pei t, 
now considered as a function of T. This is given by resummations of weak coupling 
perturbation theory at nonzero temperature. 7 )' 8 ) While all resummations fail at tem- 
peratures below Tpert j in the present view this is just because of the non-ideal terms 
in the pressure. This suggests that for all T > T c , perturbative resummations con- 
tribute only to /p ert (T). Analysis shows that / P ert(T P ert) is about 90% of the ideal 
gas value. 7 )' 8 ) If applied just to f peTt (T), it seems very possible that perturbative 
resummation might work all of the way down to the critical temperature. 7 ) With 
present day techniques, lattice simulations could test this in a precise manner in the 
pure glue theory. 

This suggests a heuristic analogy, to the operator product expansion for two 
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gauge invariant operators at short distances. Free field theory dominates as the 
distance x — > 0, like some power of x. Perturbative corrections enter with the 
same power, times a series in l/log(x), etc. Non-perturbative effects involve the 
expectation values of new operators, multiplied by the appropriate powers of x 2 . For 
the pressure, perturbative terms correct the ideal T 4 term, as a series in l/log(T), 
etc., while non-perturbative effects generate non-ideal terms, as a series in ~ 1/T 2 
times the ideal term. 

What about below T max ? I assume that a hadron resonance gas provides a 
reasonable approximation not just about T = 0, but all of the way to T c . This 
leaves T c < T < T max . For three colors in the pure glue theory, this is a nearly 
critical regime, dominated by a light excitation for the triplet Polyakov loop. Maybe 
even with dynamical quarks, this region is dominated by a light triplet loop and its 
interactions with quarks. 9 ) 

Why bother? For a pure glue theory, simulations relatively close to the contin- 
uum limit have been possible for some time. 2 ) With dynamical quarks, though, and 
in particular for the light quarks present in QCD, present day simulations are not 
close to the continuum limit. Thus (1-2) could be used by approximate models. For 
example, most hydrodynamic models 10 ) use a MIT bag model for the pressure. This 
approaches ideality much faster than a fuzzy bag model. Since hydrodynamics only 
requires the energy as a function of pressure, perhaps non-ideal corrections don't 
really matter that much, but this should be demonstrated by explicit computation. 

Terming (1-2) a "fuzzy" bag is not sheer whimsy. In the MIT bag model the 
surface of the bag is infinitely thin and fixed. 6 ) While the interface surely has nonzero 
width, it is difficult to know how to model this. For instance, if the surface of the 
bag were thin and floppy, then there would be many light excitations, in which the 
surface of the bag flops around, and the quarks remain essentially fixed. 11 ) There are 
no signs of such a multiplicity of states in the hadronic mass spectrum. 

The essential moral of the non-ideal terms in the pressure, (IT) and (1-2), is 
that the transition from a confined, to a nearly perturbative phase, is not abrupt, as 
in a MIT bag model, but gradual. This suggests that the boundary of the bag isn't 
thin, but thick. A thick bag is unlikely to flop around, since the entire width needs 
to participate. Thus light surface modes shouldn't be a problem. 

Having said this, it is not clear how to develop a more realistic bag, with a thick 
boundary. The thickness probably doesn't affect the mass spectrum of ordinary 
hadrons greatly. The analogy is still suggestive. In some loose sense nonzero tem- 
perature probes inverse distances in the QCD vacuum, T ~ 1/R: the perturbative 
vacuum emerges as T — > oo, or R — > 0; the confined vacuum, for T — > 0, or R — > oo. 
The width of the bag emerges over distances R ~ 1/T pert — ► 1/T C ; in physical units, 
for ~ 1/4 fm — » 1 fm. This is the right scale to probe the transition from going 
inside, to outside, the bag. 

§2. Wilson lines and their electric field 

The above begs the question: what is the origin of the non-ideal terms in the 
pressure of the deconfined phase? I next turn to a possible explanation in terms of 
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Wilson lines. ^ I omit all details, and many references, to concentrate on a broad 
and qualitative description. 

In resummations of perturbation theory, 7 )' 8 ) the gluon degrees of freedom remain 
the time-like, Aq, and space-like, A{, components of the vector potential. To proceed 
further, consider a straight, thermal Wilson line: 

rl/T 

ig / Aq{x, t) dr 
L{x) = P e Jo . (2-1) 

r : — > 1/T is the imaginary time, x is the spatial coordinate, and g the gauge 
coupling constant. If Aq is in the fundamental representation, L transforms as an 
adjoint matrix: under strictly periodic gauge transformations, U(x, 0) = U{x, 1/T) = 
U(x), L{x) -^U^{x)L{x)U(x). 

The trace of the Wilson line is the Polyakov loop, and is invariant under local 
gauge transformations. The motivation for considering an effective theory of the 
Wilson line rests upon measurements of the renormalized Polyakov loop. In a per- 
turbative regime, fluctuations in Aq should be small, and so suitably normalized, this 
expectation value should be near one. In contrast, if fluctuations in Aq are large, 
the expectation of the Polyakov loop is not near one. Lattice simulations in a SU(3) 
gauge theory, with or without quarks, show the (renormalized, triplet) loop is near 
one for T > T pert , and less than one for T c < T pert . 

One's first guess might be that the theory is driven into a regime of strong 
coupling. Consider, however, the "Helsinki" program of resummation. 8 -* Originally 
proposed by Braaten and Nieto, computations to four loop order were done by 
Kajantie, Laine, Rummukainen, and Schroder. The final steps are being completed 
by Laine, Schroder et a/. 8 ) They find that that the effective coupling runs with a 
mass scale ~ 2irT, so that even down to T c ~ 175 MeV, the QCD coupling is 
af (1.6 GeV) ~ 0.28. This value is not that large. 8 ) 

This suggests the perturbative construction of an an effective theory in three 
dimensions, valid over distances > 1/T. Since the renormalized Polyakov loop is not 
near unity, Aq is replaced by Wilson lines, coupled as always to the Ai. 

When I first suggested this there were several problems. 9 ) The first is how one to 
match the effective theory, at large Aq, to quantities which are computable pertur- 
batively. In the Helsinki program, this is done by computing the positions of poles 
in propagators, etc. This is fine at small Aq, but doesn't probe large Aq. 

Interfaces 12 ^ can be used to probe large Aq. These are most familiar in a pure 
gauge theory, such as SU(N). A SU(N) gauge group has a global center symmetry of 
Z(N), so that in the deconfined phase, there are N degenerate vacua: the usual per- 
turbative vacua, L = ljv, and Z(N) transforms thereof, such as L = exp(2ni/N) 1^. 

An interface interpolates between these degenerate vacua. One takes a box 
which is long in one spatial direction, say z : — > zt. At one end of the box, 
one takes one vacua, L(0) = ljv; at the other end, an inequivalent Z(N) state, 
L(zf) = exp(27ri/iV) ljy. These boundary conditions force the formation of an inter- 
face along z, which tunnels between the degenerate vacua. While the ends of the box 
represent perturbative vacua, in between one probes large Aq ~ T/g, as illustrated 
by (2-3). The amplitude for tunneling can be computed semiclassically, and gives 
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an interface tension ~ T 2 / \f~g* ■ The width of the interface is proportional to the 
inverse Debye mass, ~ l/{\f~g*T), so that a derivative expansion can be used. To 
date, computations have been carried out to ~ (\/g^) 3 ^ 2 times the result at leading 
order by Giovannangeli and Korthals Altes. 12 ' 

Such interfaces appear to be special to theories with a center symmetry, and 
so useless for QCD, where it is violated by the presence of dynamical quarks. In 
this case, however, there are other interfaces which can be used. Consider the gauge 
transformation 

Uc{T) = e 2mrTt N /N t N = ^ -(N-l) ) ' ™ 

This is only periodic up to a Z(N) rotation, U c (l/T) = exp(27ri/iV) U c (0), which is 
allowed in the absence of quarks. If one acts with U c on one end of the box, but not 
the other, then a Z(N) interface forms, because this isn't a pure gauge transformation 
in between the two ends. 

Now instead of (2-2), consider the N th power thereof, ■ This is a strictly 
periodic gauge transformation, and thus is allowed, independent of whatever matter 
fields may be present. If we act with on one end of a box, but not the other, 
what I term a U(l) interface is generated: while (L) = ljv at both ends, in between 
Aq winds around in a topologically nontrival fashion. 

To compute the properties of an interface, one needs the effective potential for 
constant Aq, which is first generated at one loop order. For a effective theory of 
Wilson lines, it is obvious to turn an effective potential of Aq into one for L. What 
stymied me 9 ' is what one does at zeroth order: how does one write the effective 
electric field in terms of Wilson lines? 

The crucial clue was provided by what appeared to be an abstruse computation. 
At one and two loop order in a SU(N) gauge theory, explicit calculation shows that 
the effective potential, computed in the presence of a large, background field for 
constant Aq — and Ai = — is invariant under the Z(N) center symmetry. 12 ' 
This is unremarkable: there is no anomaly to prevent the quantum theory from 
respecting the center symmetry of the classical theory. Diakonov and Oswald 13 ' 
then computed in the presence of background fields in which both Aq and Ai were 
nonzero, allowing Aq to be large. Assuming that the effective electric field is D^Aq, 
they found that the center symmetry appeared to be violated at one loop order. 
An error in computation seems unlikely, given that their results agree with those of 
Megias, Ruiz, and Salcedo, 14 ' who computed in the limit of small Aq and A\. 

I suggest that the problem lies in assuming that the effective electric field is 
D{Aq, and arises even at tree level. Under the large gauge transformation of (2-2), 

A dia 9 ^ A dia 9 + ^T tN A^ulMUc. (2-3) 

Hence diagonal elements of Aq are shifted by a large but constant amount, ~ Ttjy/g, 
while off-diagonal elements of Aq and Ai undergo time dependent rotations. 

In four dimensions, the original electric field is D^Aq — OqAi. The first term, 
D^Aq = diAo — ig[Ai,Ao\, changes if [Aj,tjv] 7^ 0, which always happens for some 
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(off-diagonal) components of A^. It is easy to see, however, that the time dependent 
rotation of Ai in the second term, — 8qAi, generates a similar commutator, —[Ai,t^\, 
which exactly cancels against the first term. This just reflects the fact that the orig- 
inal electric field transforms homogeneously under arbitrary gauge transformations. 

In the effective theory, the simplest guess for the effective electric field is to drop 
time derivatives, and take it to be D{Aq. While valid at small j4o, 7 )' 8 ) this can't be 
right at large Aq. The above argument shows that D{Aq is not invariant under the 
large gauge transformations, IA C , which enforce the center symmetry. 13 ) Similarly, 
if we take , the effective theory is not even invariant under large, but strictly 
periodic, gauge transformations. 

The resolution is to construct the effective electric field from Wilson lines. These 
transform like a phase under IA C: and are invariant under . To leading order in 
weak coupling, the effective electric field is 

Ef(x) = — LUx) Di L(x) . (2-4) 

This is shown by demonstrating that the interface tensions, either Z(N) or U(l), agree 
between the effective and original theories. Eq. (2-4) does change beyond leading 
order, and is multiplied by other gauge invariant terms, such as |trL| 2 , |trL 2 | 2 , etc.; 
with coefficients which begin at ~ g 2 . 

It is well known that the mapping between the fields in an effective theory, and 
the original theory, is indirect. Usually, however, one only sees this at next to leading 
order in some expansion, so the differences are small. In the present instance, because 
one is constructing an effective theory for large Aq, it arises even at zeroth order. 
The problems found before 13 ) are presumably solved by matching to an effective 
theory constructed from the effective electric field, and not functions of D^Aq. 

To leading order, the effective Lagrangian is that of a gauged, principal chiral 
field: 

^L.^.. L) = i tr G% + ^ tr \L^D t L\ 2 . (2-5) 

This is nonrenormalizable, but this shouldn't be a problem, since the effective the- 
ory is only valid over distances > 1/T. A related linear model, constructed to be 
renormalizable, has been analyzed by Vuorinen and Yaffe. 15 ) 

On the lattice, the Lagrangian of (2-5) was first written down by Banks and 
Ukawa, 16 ) as the simplest kinetic term for an adjoint field, L. It is not obvious that 
(2-5) applies in the continuum, even at leading order in weak coupling. It should 
be possible to construct effective Lagrangians to the same order as interfaces, in the 
original theory, have been computed. 12 ) 

One can show that the four dimensional instanton number equals the winding 
number of the Wilson line. 1 ) Dai and Nair 17 ) showed that non-abelian hydrodynamics 
has color Skyrmions. This suggests that the effective model might have electric 
Skyrmions: 17 ) solutions stabilized by a nonzero winding number for L, and yet which 
are not instantons, since they only carry electric, and not magnetic, fields. 

How does the deconfining transition arise? For the potential for L computed 
perturbatively, order by order the deconfined vacuum is always proportional to the 
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unit matrix, (L) ~ Ijv- This could be possible all of the way to T c — but then one 
would not expect non-ideal terms in the pressure. To destabilize the perturbative 
vacuum, it is necessary to add, by hand, a non-perturbative mass term for the 
Wilson line. The simplest example is ~ T 2 Bf |trL| 2 . The notation Bf, and the 
T 2 , is motivated by the previous section; adding such a mass term is standard in a 
Landau-Ginzburg type of analysis. 

With such a term, (L) ^ 1. In Ref. 1, I argued that at infinite N c , the confined 
vacuum is characterized by complete eigenvalue repulsion; also, that lattice results 
suggest that this is approximately true at small N c . The appearance of eigenvalue 
repulsion is clear when the spatial volume is small, as when the theory lives on a very 
small sphere. At infinite N c , Aharony et. al. ls ^ showed that the effective theory only 
involves a constant mode, and reduces to a random matrix model for that mode. As 
typical of random matrix models, eigenvalue repulsion arises from the Vandermonde 
determinant in the measure of the group integral. 

In infinite volume, instead of one random matrix, there is a field theory of (not 
quite) random matrices. Such field theories have been studied little. The inter- 
esting, and gauge invariant, question is how confinement arises from the dynamical 
generation of eigenvalue repulsion. 
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